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ABSTRACT 

We have undertaken a search for the infrared emission from the intracluster dust in the Coma 
cluster of galaxies by the Multiband Imaging Photometer for Spitzer. Our observations yield the 
deepest mid and far-infrared images of a galaxy cluster ever achieved. In each of the three bands, 
we have not detected a signature of the central excess component in contrast to the previous report 
on the detection by Infrared Space Observatory (ISO). We still find that the brightness ratio between 
70/zm and 160/im shows a marginal sign of the central excess, in qualitative agreement with the ISO 
result. Our analysis suggests that the excess ratio is more likely due to faint infrared sources lying 
on fluctuating cirrus foreground. Our observations yield the 2a upper limits on the excess emission 
within 100 kpc of the cluster center as 5 x 10~ 3 MJy/sr, 6 x 10~ 2 MJy/sr, and 7 x 10~ 2 MJy/sr, 
at 24, 70, and 160 /im, respectively. These values are in agreement with those found in other galaxy 
clusters and suggest that dust is deficient near the cluster center by more than 3 orders of magnitude 
compared to the interstellar medium. 

Subject headings: galaxies: clusters: general - galaxies: clusters: individual (Abell 1656) - infrared: 
general - intergalactic medium 



1. INTRODUCTION 

The presence of dust grains in the diffuse intergalac- 
tic medium is still an open question. If they exist, they 
should have great impacts not only on our understanding 
of galaxy evolution but also on interpretation of high- 
redshift observational data. From X-ray measurements 
of heavy metal lines, it is evident that clusters and groups 
of galaxies contain a sig nificant amount of ele ments pro- 
cessed in galaxies (e.g., iTamura et al.l I2004L and refer- 
ences therein). Such elements are likely to have been 
expelled out of galaxies through ram pressure stripping, 
galaxy mergers and interactions, blowout by supernovac, 
etc. Dust may also be expelled by these processes or 
by some i ndependent mechanism such as radiation pres- 
sure (e.g., iChiao fe Wickramasinghe1ll972l : lFerrara et alJ 
Il991f> . Th e X-ray emit ting hot gas can collisionally heat 
the dust (|Dwek |[i986f ). while it can also destroy small 
dust grains via sputte ring (|Draine fc Salpeter l fl979l : 
ITsai fc Mathewsl fl995h . Dust grains coexistent with 
the X-ray emitting gas have indee d been detected in 
a number of elliptical galaxies (e.g., iKnapp et al.l 119891 : 
iKaneda et alll2007t ITemi et al.ll2007[ ). The temperature, 
the amount, and the size distribution of the dust in hot 
plasma can be quite different from those in the normal 
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interstellar medium (see Dwek & Arendt 1992 for re- 
view on dust-gas interactions). Once detected firmly, 
the intracluster or intragroup dust will definitely provide 
a powerful probe of galaxy evolution and the dust-gas 
interaction. 

There have been several suggestions and debates re- 
garding the presence of dust i n cluster s and groups of 
galaxies. Over five decades ago. lZwickv I (|1957[ ) first sug- 
gested an enhanced extinction to ward the Coma clus - 
ter of galaxies. More recently, iGirardi et al.l (|1992l ) 
pointed out that the redshift asymmetries of mem- 
ber galaxies in nearby groups are consistent with the 
presence of dust in the intragroup medium. The ob- 
served oxygen K edge in an X-ray spectrum of the 
Perseus clust er can also be attribu t ed to intracluster 
dust grains (| Arnaud fe Mushotzkvl 119981 ). The ex- 
tended submiUimeter emission detected in a rich galaxy 
clust e r may partly be due to dust (jKomatsu et al.l 
119991: IKitavama et alJ I2004D . There have been ex- 
tensive searches for the enhanced visual ext inction 
towa r d high redshift objects behind clusters ( | Maoz I 
19951: INollenberg et all 120031: IMuller et all |2008|) and 
Chelouche et al.l (|2007[ ) reported the detection of redden- 
ing toward z ~ 0.2 clusters. On the other hand, the in- 
traclust er dust is like ly to enhance the efficiency of gas 
cooling (|Dwek 1 119871 ) that is in apparent contradict ion 
with recent X-ray observations (Peterson et al. 20fJ3l). 

More direct evid ence has been searc hed for in the 
far- infrared bands. Hickson ct al. (1989) reported that 
the far-infrared emission observed with IRAS is en- 
hanced by about a factor of 2 in compact groups of 
galaxies compared with a sample of isolated galax- 
ies. The extended excess emissio n is also inferre d 
toward Abell 262 and Abell 2670 (jWise et al.l [l99l . 
ISulentic fc De Mello R abacal (119931) .howev er, pointed 
out that the results of IHickson et alJ (|1989f ) are likely 
to have been overestimated mai nly due to t h e lim - 
ited spatial reso l ution of IRAS. IStickel et all (|1998f ) 
and IStickel et alJ (|2002fl reported the detection of ex- 
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tended excess far-infrared emission from the Coma clus- 
ter, which can be attributed to thermal emission from 
the intracluster dust. They find such emission only in 
one out of six cluste rs they observed w i th Inf rared Space 
Observatory (ISO). iMontier fe Giard I ([2005) combined 
the IRAS maps toward 11,507 galaxy clusters and statis- 
tically detected the excess emission, although the contri- 
bution of member galaxies is yet to be fully removed. 

Major difficulties in the previous infrared searches are 
the limited sensitivity and spatial resolution of the detec- 
tors. These limitations may have led to significant con- 
tamination of individual galaxies as well as Galactic cir- 
rus confusion. It will therefore be meaningful to perform 
infrared observations using current facilities with much 
impro ved sensitivity and spatial resolution. iBai et aTl 
(2007) searched for the intracluster dust emission from 
Abell 2029 with Spitzer and reported that the cirrus noise 
limits the detection of the intracluster dust component. 

In this paper, we present the results of our deep ob- 
servations toward the central region of the Coma cluster 
(Abell 1656) with the Multiband Imaging Photometer 
for Spitzer (MIPS) on board the Spitzer Space Telescope. 
Coma lies at high galactic latitude (b — 88°) where the 
cirrus level is expec ted t o be among the low est. In fact, 
iStickel et all {1998) and lStickel et all (|2002t ) found with 
ISO the excess far-infrared emission toward Coma based 
on the brightness ratio /120/iia/ Ixsoumi while the S/N at 
each wavelength is still insufficient to detect the excess di- 
rectly and the contamination of galaxies is not removed. 
We attempt to probe the excess emission more directly by 
approaching the confusion limits of extragalactic sources 
and Galactic cirrus in each band. 

Throughout the paper, we assume a standard set of 
cosmological parameters, O m = 0.27, Qa = 0.73, and 
h = 0.71 consistent with the WMAP five year data 
(jKomatsu et al.ll2009f h An angular size of 1" corresponds 
to a physical size of 0.46 kpc at the redshift z — 0.0231 
of the Coma cluster. 

2. OBSERVATIONS AND DATA REDUCTION 

The central region of the Coma cluster was observed at 
24, 70 and 160 /im with the MIPS in our G03 program 
on 2007 January 19. Total duration of the observation 
was 41 ksec. Linear scans were performed in a slow scan 
mode covering the 5' x 50' region across the cluster as 
shown in Figure [T] (the scan width was reduced to ~ 
2.5' at 70 /im due to dead pixels). The scan paths were 
carefully chosen so as to avoid bright galaxies and to 
lie close to the cluster center; the scan center was fixed 
at the mid-point (12 h 59 m 51.9 s , +27°58'05") between two 
massive central galaxies, NGC4874 and NGC4889 with 
the position angles of -25.4°, -28.9° and -25.2° at 24, 
70 and 160 /jm, respectively. The above scan center is 
3.2' northeast of the X-ray center of the cluster in the 
ROSAT PSPC image (12 h 59 m 40.1 s , +27°56'56"). This 
separation is small compared to the extent of the X-ray 
emission, which is characterized by the core radius of 10'. 
In total, 32 scans were performed in each band, 16 of 
which repeatedly covered a 5' x 30' region on cither side 
of the scan center. The integration time per pixel after 16 
scans was 1600, 1600 and 160 sec at 24, 70 and 160 /im, 
respectively. Two sets of the 16 scan legs overlapped and 
the integration time was doubled in the central 5' x 10' 
region, where the X-ray emission is the most prominent. 



The data were first processed with the MIPS pipeline, 
version S17 (|Gordon et al.ll2005f) to yield the Basic Cali- 
brated Data (BCD). To further improve the images, the 
pixel-dependent and scan-dependent backgrounds were 
corrected as follows. In doing so, special care was taken 
to maintain the extended nonlinear component of the 
emission. All the results regarding the extended emis- 
sion are based on unfiltered images. 

First, we create a median fiat field of detector pix- 
els by stacking the normalized BCD frames, excluding 
the regions within 10' from the cluster center. The nor- 
malization is done by dividing each BCD by its median. 
This is intended to separate any large-scale variations, 
e.g., due to the zodiacal light, from the pixel-dependent 
backgrounds and applied solely to create the flat field. 
We have checked that the flat field so obtained is prop- 
erly normalized; its mean is equal to unity within 0.2%. 
We then divide the original BCD (not normalized) by 
the normalized median flat field. At 24 /im, we further 
remove scattered light that produces a background mod- 
ulation as a function of the scan mirror position; the 
average offset at each scan mirror position is computed 
excluding the data within 10' from the cluster center. As 
mentioned below, these processes will be re-run after the 
source identification. 

Second, we create for each band a mosaic im- 
age using the MOPEX software package, version 18.2 
(jMakovoz fc Khan"! I2005f ) developed at the Spitzer Sci- 
ence Center. The pixel size is taken to be a nominal 
value of 2.45", 4.0" and 8.0" at 24, 70 and 160 fim, re- 
spectively. We exclude the first five BCD frames in each 
scan leg to omit boost frame transients. We also ex- 
clude the frames affected by stimulator flashes at 70 and 
160 /im. The background matching is applied after com- 
bining the frames in each scan leg; we only correct the 
constant offsets among 32 scan legs to assure sufficiently 
large overlapping area and to preserve any extended fea- 
ture of the emission. 

Third, we identify point sources in the mosaic images 
down to the la detection limit (~ 0.02, ~ 2 and ~ 20 
mJy at 24, 70 and 160 m respectively) using APEX 
(jMakovoz fc Marleau 2005). This limit corr esponds to 
the 2 ~ 3cr confusion level in each band (|Dole et al.l 
2004b). The regions around the sources are masked out 
and unused in our analysis. More specifically, the source 
identification and the masking procedure are done on 
the basis of median-subtracted images; a median-filtered 
image is created and subtracted from the unfiltered im- 
age solely for this purpose. To improve the efficiency 
of source detection against unknown background fluc- 
tuations, we vary the size of a median filter box from 
5 x 5 to 25 x 25 pixels. For every source identified in 
each median-subtracted image, we estimate the radius 
T'mask at which the surface brightness, including that of 
the side-lobes, drops below a certain threshold noise level 
using the Point Response Function. We adopt 0.5cr as the 
threshold noise level at 24 and 70/im, and la at 160 /im, 
considering the lower uniformity of the data in the latter 
band. We then return to the unfiltered image and mask 
all the pixels within r mas k from the sources identified in 
any of the median-subtracted images. We have checked 
that the correlation coefficient between the mean bright- 
ness of masked pixels (r < r mas k) and that of unmasked 
pixels at r mas k < r < 2r mas k is less than 0.2 in any of 
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the median-subtracted images. We thus consider that 
our masking procedure successfully removes the contri- 
bution of resolved sources. 

Strictly speaking, the source identification and the 
background correction are not independent of each other. 
We therefore repeat the above procedure from the flat- 
field correction excluding the pixels around the point 
sources in the BCD. 

After applying the above corrections, we subtract from 
each unfiltered image the linear gradient over the entire 
map, so as to remove the variation on scales larger than 
the galaxy cluster. Such a gradient is greatest at 24 
/im due to the zodiacal light; it causes the difference of 
A/„ > 1 MJy/sr between the map edges which is more 
than 2 orders of magnitude larger than the average noise 
level of the image (Table [TJ). By simply subtracting the 
linear component, we avoid relying on specific models 
of the zodiacal light emission that are still uncertain at 
scales probed here. To preserve the average brightness of 
the entire map, the brightness at the scan center is kept 
unchanged by the subtraction. 

3. RESULTS 

Figure [2] shows the mosaic images after the subtraction 
of the linear gradient. The point sources are explicitly 
shown for display purposes. The pixel sizes of the im- 
ages, the FWHMs of the Point Response Function, and 
the average noise levels are listed in Table [TJ We have 
assigned a noise to an individual pixel as the larger of 
the uncertainty estimated in the MIPS pipeline and the 
standard deviation of actual measurements at that pixel. 
The latter is on average twice as large as the former at 
70 and 160 /im, and both are comparable at 24 /im. This 
choice should enable us to take account of any residual 
systcmatics in each pixel more properly. Note that the 
sensitivity at 160 /im is likely limited by Galactic cirrus; 
the lcr fluctuation level expected at the position of Coma 
is 5. 4 x 10~ 3 MJy/sr at 70 /im and 0.20 MJy/sr at 160 
/xm (peong et al.l(2005D . 

Masking the regions around the detected sources re- 
moves 38%, 66% and 64% of the pixels in the 24, 70 and 
160 /im images, respectively. We plot the histograms 
of the pixels with and without masking in Figure [3J In 
all bands, the histograms before masking are positively 
skewed due to bright sources. Once the sources are ex- 
cluded, the data follow the Gaussian distribution more 
closely, while there are small excesses at high a. These 
excesses are due to low-significance data affected by noisy 
detector arrays or low coverage. Since the noise assigned 
to such data is large, their contribution to the mean 
brightness is negligible. At 160 /im, the la range of the 
distribution after the source mask, 0.19 MJy/sr, is likely 
affected by the flu ctuations of Galactic cirrus emission 
()Jeong et al.ll2005f ). At 24 and 70 /im, the la ranges, 
0.011 and 0.11 MJy, are comparable to the average noise 
level listed in Table [TJ 

Figures [4][6] show the surface brightness distribution 
along the scan path as a function of the distance from 
the X-ray center of the cluster (negative values of the dis- 
tance are assigned to the pixels at lower declination than 
the X-ray center). We plot the brightness at each pixel 
and their binned average. Unless otherwise stated, the 
error bar indicates the la statistical error of the mean in 
each bin. The 70 /im data are lacking within 1.7' from 



the X-ray center due to dead pixels. We hence adopt 
fiducially the bin size of 3.6' (100 kpc) in the present pa- 
per so that it fully covers the deficit region while it is still 
sufficiently smaller than the X-ray core radius of 10'. In 
all bands, the fluctuations in the brightness, particularly 
the spikes due to point sources, are significantly reduced 
after the source mask, yet there is no indication of the 
central excess. 

In the bottom panels of Figures [4]-[6j we plot the mean 
brightness of the removed component, i.e., the difference 
between the unmasked and masked brightness. From 
our linear scans within ~ 30' from the cluster center, 
we do not find a central concentration of the integrated 
source intensities either. The average levels of the re- 
moved brightness are 0.012 MJy/sr, 0.093 MJy/sr and 
0.16 MJy/sr at 24, 70 and 160 /im, respectively. 

We further combine the brightness in the two directions 
from the cluster center and compare its radial profile with 
model predictions in Figure [7J The prediction for the in- 
tr acluster dust emission i s based on the theoretical model 
of lYamada fc Kitavama I (|2005h (see the next section for 
details). For direct comparison with the ISO results, we 
have adjusted the total amount of dust to reproduce the 
excess brightness of 0.2 MJy/sr at 120 /im within the 
radius of 1.5' inferred in iStickel et~ai1 ([2002). bmce we 
are only interested in the difference of the signal across 
the cluster, an offset is added to the model brightness 
to match the mean brightness of the whole map at the 
edge. At 24 /im, the absence of the excess is as expected 
since the predicted temperature of the intracluster dust 
is relatively low, 20-30K (see also Dwek et al. 1990). At 
70 and 160 /im, the foreground/background fluctuations 
are comparable to the expected levels of the emission and 
we do not find a sign of the central excess in either band. 

From the observed brightness within r = 3.6' relative 
to the mean of the entire map, we obtain 2a upper lim- 
its on the central excess as 5.2 x 10~ 3 , 6.0 x 10 -2 and 
6.6 x 10 -2 MJy/sr, at 24, 70 and 160 /im, respectively. 
In deriving these limits, we have taken into account the 
background/foreground fluctuations that dominate over 
the statistical noise as follows. We have randomly se- 
lected a circular region of radius 3.6' from the map in 
each band 100 times. From the variance of the mean 
brightness of selected regions, we have assigned an la- 
error of 1.9 x 10~ 3 , 2.1 x 10~ 2 and 6.6 x 10~ 2 MJy/sr at 
2 4, 70 and 160 / im, re specti vely. 

IStickel et all (|1998l ) and IStickel et all (j2002fl found 
with ISO the excess far-infrared emission toward the cen- 
ter of this cluster based on the surface brightness ratio 
^i20Aim/-^i80Mm- Although the scan paths and the ob- 
served bands are not exactly the same in our observa- 
tions, we attempt to mimic their analysis by presenting 
Izifna/Iroum and hoia*/ hwpan in Figure M While the 
fluctuation is large, there appears to be a weak sign of a 
central enhancement in Ijonm/Ii60tJ,ni- The enhancement 
is largely due to a variation in the 160 /im map around 
the cluster center (Figj2|); there are bright extended re- 
gions at ~ 10' on each side, presumably produced by far- 
infrared sources lying on fluctuating cirrus foreground. 
We have checked that no counterparts are detected in 
the Chandra X-ray image of this region. While some 
of the sources are apparent and removed at 70 /im, the 
lower spatial resolution and the larger cirrus fluctuation 
prevent us from removing them completely at 160 /im. 
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In turn, the cluster center appears to be darker than its 
surrounding at 160 /im, resulting in the apparent excess 
in /70/im/^i60/im- Figure [6] also indicates that the depres- 
sion at 160 /im is not due to over-removal of the sources 
at the center. We therefore do not consider this to be a 
detection of the intracluster dust. 

4. DISCUSSION 

While there is no direct evidence of the central excess in 
any of the three bands observed by the MIPS, we still find 
a weak sign of excess in the brightness ratio IjQ^m/Iimimx 
th at may resemb le the exces s in Ji20um/^i 8 pmii r eported 
by IStickel et ail (jl998fl and IStickel et ail (|2002f ). It is 
therefore meaningful to explore the consistency of these 
results. 

For this purpose, we subtract from our data the aver- 
age zodiacal light and the cosmic infrared background 
levels estimat ed by the Spitzer background estimator 
(jReach 112001 ) at the scan center position for the date of 
observation; J ZO d = 7.21 and 1.41 MJy/sr, J C ib = 0.21 
and 1.30 MJy/sr at 70 and 160 /im, respectively. We 
subtract the above values from the whole map since we 
have already corrected the linear gradient of the map 
relative to the scan center. We have checked that sub- 
tracting position-dependent zodiacal light levels prior to 
the correction of the linear gradient changes the result- 
ing brightness ratio by no more than 5%. The uncer- 
tainties in the absolute calibration are not significant ei- 
ther; 5% at 70 nm and 12% at 160 /im (|Gordon et all 
l2007t IStansbe rrv et al.ll2007[ ). However, simply because 
the mid-infrared brightness is dominated by the zodiacal 
light, the nominal error of ~ 15% in its model (jReach I 
2000) will change the brightness after subtraction at 70 
and 160/im, - 1.1 MJy/sr and ~ 1.9 MJy/sr, by ~ 100% 
and ~ 10%, respectively. Keeping these large uncertain- 
ties in mind, we investigate the origin of the excess in 
the band ratio as follows. 

Figure [9] shows the radial distribution of the bright- 
ness ratio ^7o" m /^i60Aim after the correction mentioned 
above. The error bars include only the statistical uncer- 
tainties. Assuming that the corrected brightness distri- 
bution is produced by a combination of a central excess 
and more extended component, e . g., Ga lactic cirrus, we 
apply equation (7) of lStickel et al.1 (|2002t ) to our observed 
bands and obtain I^ c ^ s ~ 6.4 x 10~ 2 MJy/sr within the 
radius of 3.6'. This value is comparable to the 2a upper 
limit of 6.0 x 10 -2 MJy/sr derived solely from the 70 /im 
data in the previous section. The central excess, if real, 
should have been nearly canceled at 70 /im by the fluctu- 
ation of the other component. The interpretation of the 
central excess for the origin of the observed Irofim/Iieo^m 
profile is therefore rather unlikely, although it is not fully 
rejected. 

Note that equation (7) of lStickel et alj (|2002D is based 
on a key assumption that the excess component con- 
tributes only to the variation of the shorter wavelength 
data. Alternatively, it is also possible to interpret the 
same data with the variation of the longer wavelength 
data. In the present observation, a plausible explanation 
may be that the corrected brightness is still dominated 
by unresolved sources lying on the fluctuating cirrus fore- 
ground, i.e., the absolute values of the signal in both 
bands are dominated by such sources while the fluctua- 



tions at 160 /im are mainly produced by underlying cirrus 
emission. In fact, the value of Iffi/Im/ Il&aum = 0-55 ~ 0.6 
shown in Figure [9] corresponds to the dust temperature 
of Td ~ 25 K for the absorption efficiency of Q a bs oc A~ 2 
(Drainc & Lee 1984). This temperature is higher than 
the nominal temperature of Galactic cirrus T4 ~ 18 K 
for which the brightness ratio is lower by more than a 
factor of 5. The mean values of IfoJf m ~ 1.1 MJy/sr and 
^I60uin ~ 1-9 MJy/sr are also higher than the levels of 
the cirrus emission expected for this region, suggesting 
the contribution of yet unresolved components. 

From the limits on the infrared brightness in each 
band, we further attempt to derive constraints on the 
amount of the int racluster dust. Model pre dictions 
shown in Figure [7] (Yamad a k. Kitavama ll2005f ) rely on 
an assumption that the dust is supplied continuously to 
the intracluster medium from galaxies with the initial 
size distribution sim ilar to our Galaxy, dn lmt /da oc a~ 3,5 
(jMathis et al.1 11977D . and has reached a steady state 
against sputtering with the size distribution dn/da oc 
a -2.5 g a j ax y distribution is assumed to trace that 

of the dark mat ter given by the Nav arro-Frenk -White 
density profile (|Navarro et al.l fl997). and the gas den- 
sity follows the isothermal beta model. Given these, the 
curves plotted in Figure [JJ correspond to the dust sur- 
face density of S</ = 1.1 x 10 3 M Q kpc -2 within the 
central 3.6' of the cluster. The 2a upper limits on the 
excess brightness derived in the previous section are 1.3 
and 0.37 times the model brightness averaged within the 
central 3.6' at 70/im and 160/iin, respectively. Consider- 
ing much lower uniformity of the data at 160 /im, we take 
the limit of the 70/im data to obtain E<j < 1.4 X 10 3 M 
kpc -2 . This limit corresponds to the dust-to-gas-mass 
ratio of less than 1 x 10~ 5 within the central 100 kpc, 
which is nearly 3 orders of magnitude smaller than the 
Galactic value (jDraine 1120031 ) . 

The average visu al extinction A y is related to the dust 
surface density as (|Spitzer 1119781 ) 

Ay = 0.011 ( s ) 

v 2 j V10 3 M Q kpc- 2 y 
x( pd VV-g-V 1 (i) 

V3 g cm- 3 / V0-1/W 

where Q cx t is the extinction efficiency factor, pd is the 
dust mass density, and a is the grain radius. Adopt- 
ing the standard values for Qext, pd and a quoted in 
the above equation, our derived limit of < 1.4 x 10 3 
M Q kpc -2 corresponds to A y < 0.016 within the cen tral 
100 kpc. On the other hand, iChelouche et al.l (|2007| ) re- 
port the detection of the mean E(B — V) of a few times 
10~ 3 mag toward S PSS QSOs behind ~ 10 4 clusters, and 
IMuller et all (|2008l) obtain the limit (Ay) < 0.024 (2a) 
from reddening of ~ 90, 000 background galaxies for 458 
clusters. While these optical measurements are given for 
much larger (> Mpc) scales, our result is still consistent 
with them. 

Finally, gravitational lensing can impact on the emis- 
sion of background sources behind galaxy clusters. As 
the total surface brightness is conserved by lensing, the 
bias in the resolved or unresolved brightness is caused 
solely by the sources tha t exceed the detection li mit as 
a result of magnification (|Refregier fc Loebl [l997'). The 
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number of such sources is / = //* — 1 times the number 
of unlensed sources already lying above the flux limit 
S, where \x is the magnification and a is the index of 
the unlensed source count, N(> S) oc S~ a . The frac- 
tional change in the resolved brightness has the same 
sign as and is less than /, because the sources con- 
tributing to the change have the lowest flux among those 
resolved. Note that the factor / is different from the 
bias in the overall surface dens ity given by ii Q_1 (e.g., 
IBartelmann fc Schneider! 12001). Recent surveys with 
Spitzer yi eld the mean values o f a = 0.5 down to 35 /xJy 
at 24 fim (IPapovich et al.ll2004D an d a = 0.6 down to 1.2 
mJy at 70 /xm (|Fraver et al.ll2006l ). The depths of these 
measurements are comparable to those of the present ob- 
servation. At 160 /im, the published counts are currently 
available for S > 50 mJy, and an extrapolation using the 
mode l oflLagache et alt (12004ft yields a ~ 1.5 at S ~ 20 
mJy ([Dole et al.l|20 04a). In the case of Coma, the Ein- 
stein radius is ~ 30" and the magnification is expected 
to vary from fi ~ 1.2 at r = 3.6' to n ~ 1.02 at r = 30' 
for a singular isothermal mass profile. Taken together, 
the above factor in this radius range varies from / ~ 0.1 
to 0.01 at 24 fim and 70 fim, and / ~ 0.3 to 0.03 at 160 
/im. The fractional change in the brightness of the re- 
moved component shown in Figures 4-6 should therefore 
be smaller than the variations along the scan path. There 
is in fact no sign of systematic increase in the resolved 
brightness toward the cluster center in any band. Owing 
to the conservation of the total brightness, the bias in the 
residual unresolved background should also be small. We 
hence consider that the impact of gravitational lensing is 
negligible in the present results. 

5. CONCLUSIONS 

In this paper, we have undertaken a search for the in- 
frared emission from the intracluster dust in the Coma 
cluster of galaxies by the MIPS on board the Spitzer 
Space Telescope. Our observations yield the deepest mid 
and far-infrared images of a galaxy cluster ever achieved 
with average ltr sensitivities per pixel of 1.3 x 10 ~ 2 
MJy/sr, 1.1 x lO -1 MJy/sr and 1.3 x 10" 1 MJy/sr at 
24, 70 and 160 /xm, respectively. In each of the three 



bands, we have not detected a signature of the central 
excess component. Our observations yield the 2a up- 
per limits on the excess emission within 3.6' (100 kpc) 
of the cluster center as 5.2 x 10 -3 MJy/sr, 6.0 x 10~ 2 
MJy/sr, and 6.6 x 10" 2 MJy/sr, at 24, 70, and 160 /jm, 
respectively. These results ar e consistent with the values 
reported bv lBai et all (j2007l ) for Abell 2029, as well as 
those derived from a sta cking analysis of 11,507 clusters 
IjMontier fc Giard ll27)05ft . 

With the improved spatial resolution and sensitivity 
of the MIPS, we have obtained much more robust limits 
on the excess emission from this cluster than previous 
observations. We still find that the brightness ratio be- 
tween 70/xm and 160/xm shows a marginal sign of the 
ce ntral excess, in qual itativ e agreement with th e finding 
of lStickel etall (fl998ft and iStickel elaTl ([20021 ) by ISO. 
Although we cannot fully reject the contribution of the 
intracluster dust, it is more likely due to unresolved in- 
frared sources lying on fluctuating cirrus foreground. 

Combined with t he theoretical model of 
lYamada fc Kitavama I (|2005l ). we have derived a 
limit on the surface mass density of the intracluster dust 
as < 1.4 x 10 3 M kpc -2 . This limit corresponds 
to the dust-to-gas-mass ratio of less than 10 -5 within 
the central 100 kpc, which is nearly three orders of 
magnitude smaller than the Galactic value, implying 
that the dust destruction is very efficient at the center 
of the cluster. Our result can further be converted to 
the limit on the visual extinction of Ay < 0.02, which 
is consistent with the results derived from statistical 
reddening of backg round objects behind a large sample 
of ga laxy clusters (jChelouche et al.l 120071 : iMuller et al.1 
12008ft . 
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Fig. 1. — Fields-of-view of our Spitzer/MIPS observations at 24 fim (black box), 70 (im (white box) and 160 /an (green box) overlaid 
on the ROSAT PSPC X-ray image of the Coma cluster. The X-ray image is scaled logarithmically in units of photons/sec/pixel in the 
0.2—2 keV band, after smoothing by a Gaussian filter with a = 4 pixel, where the pixel size is 14.9 arcsec. The center of the X-ray surface 
brightness is indicated by a circle, and the positions of two massive central galaxies, NGC4874 and NGC4889, are marked by a cross and 
a plus, respectively. 
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Fig. 2. — Mosaic images at 24 /im (left panel), 70 fim (middle panel), and 160 (im (right panel), respectively. The positions of two massive 
central galaxies, NGC4874 and NGC4889, are marked by a cross and a plus, respectively. 
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Fig. 3. — Histograms of the image pixels before (dashed curves) and after (solid curves) the source mask at (a) 24 fim, (b) 70 fim, and 
(c) 160 fim, respectively. Dotted curves indicate, for reference, the Gaussian distribution with a = 0.011 MJy/sr, 0.11 MJy/sr, and 0.19 
MJy/sr at 24/im, 70/xm, and 160/^m, respectively. 
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Fig. 4. — Top: distribution of the 24 fim brightness along the scan path before the source mask. For clarity, we randomly select and 
display 1/10 of all the data points. Solid curve indicates the mean brightness with the bin size of 3.6' (100 kpc). Middle: same as the top 
panel except that the surface brightness after the source mask is shown. Bottom: dashed curve shows the binned average of the removed 
component, i.e., the difference between the mean brightness plotted in the top and middle panels. For comparison, solid curve is the same 
as that in the middle panel except that the average brightness of the whole map is subtracted and the vertical scale is enlarged. In all the 
panels, the distance is measured from the X-ray center of the cluster with negative values indicating the direction of decreasing declination. 
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Fig. 5. — Same as Figure [4] except that the 70 /an brightness is shown and, for clarity, 1/2 of all the data points are displayed in the top 
and middle panels. The 70 (im data are lacking within 1.7' from the X-ray center due to dead pixels. 
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Fig. 6. — Same as Figure [4] except that the 160 /im brightness is shown and all the data points arc displayed in the top and middle panels 
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Fig. 7. — Radial profiles of the surface brightness after the source mask at 24 (im (top panel), 70 /im (middle panel), and 160 fim (bottom 
panel), respectively. Solid and dashed curves indicate the mean brightne ss with the bin size of 3.6' a nd 7.2' , respectively. For reference, 
the prediction for the intracluster emission based the theoretical model of Yamada & Kitayama (2005) is plotted in dotted curves (see the 
text for details). 
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Fig. 8. — Distribution of the brightness ratios /24/jm /^70Mm (top panel) and /70fim/^l60/jm (bottom panel) after the source mask along 
the scan path. The distance is measured from the X-ray center of the cluster with negative values indicating the direction of decreasing 
declination. Solid and dashed curves indicate the brightness ratios with the bin size of 3.6' and 7.2', respectively. For clarity, the dashed 
curves have been shifted downwards by 0.1 in both panels. 




Fig. 9. — Radial profile of the brightness ratio ^7o" m /^i60(zm a fter the source mask and the subtraction of the zodiacal light and the 
cosmic infrared background. The error bars include only the statistical uncertainties. Solid and dashed curves indicate the brightness ratios 
with the bin size of 3.6' and 7.2', respectively. 



